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Abstract 

O 

Warm inflationary universe models in the context of a Chaplygin gas equation are studied. 



(N 



General conditions required for these models to be realizable are derived and discussed. By using 
a chaotic potential we develop models for a dissipation coefficient of the form r oc (ft n , with n = 
or n ^ 0. 
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I. INTRODUCTION 



It is well know that warm inflation, as opposed to the conventional cool inflation, presents 

n 

the attractive feature that it avoids the reheating period [jj. In these kind of models dissipa- 
tive effects are important during the inflationary period, so that radiation production occurs 
concurrently together with the inflationary expansion. If the radiation field is in a highly 
excited state during inflation, and this has a strong damping effect on the inflaton dynamics, 
then, it is found a strong regimen of warm inflation. Also, the dissipating effect arises from 
a friction term which describes the processes of the scalar field dissipating into a thermal 
bath via its interaction with other fields. Warm inflation shows how thermal fluctuations 
during inflation may play a dominant role in producing the initial fluctuations necessary 
for Large-Scale Structure (LSS) formation. In these kind of models the density fluctuations 
arise from thermal rather than quantum fluctuations These fluctuations have their origin 
in the hot radiation and influence the inflaton through a friction term in the equation of 
motion of the inflaton scalar field {sj]. Among the most attractive features of these models, 
warm inflation end at the epoch when the universe stops inflating and "smoothly" enters 
in a radiation dominated Big-Bang phase [1]. The matter components of the universe are 
created by the decay of either the remaining inflationary field or the dominant radiation 
field Q. 

On the other hand, the generalized Chaplygin gas has been proposed as an alternative 
model for describing the present accelerating of the universe. The generalized Chaplygin 
gas is described by an exotic equation of state of the form 

Pch = -4-, (1) 

Pch 

where p c h and p c h are the energy density and pressure of the generalized Chaplygin gas, 
respectively. (3 is a constant that lies in the range < (3 < 1, and A is a positive constant. 
The original Chaplygin gas corresponds to the case /3 = 1 jfij]. The above equation of state 
leads to a density evolution in the form 5] 



Pch 



a B 



i 



1+/3 

(2) 



a 3(l+/3) 

where a is the scale factor and B is a positive integration constant. 

The Chaplygin gas emerges as a effective fluid of a generalized d-brane in a (d+1, 1) space 



time, where the action can be written as a generalized Born-Infeld action These models 
have been extensively studied in the literature |7|. The model parameters were constrained 
using current cosmological observations, such as, CMB and supernova of type la (SNIa) 

In the model of Chaplygin inspired inflation usually the scalar field, which drives inflation, 
is the standard inflaton field, where the ener gy d ensity given by Eq. ([2|) , can be extrapolate 
for obtaining a successful inflationary period 1 Of ] . Recently, tachyon-Chaplygin inflationary 



universe model was studied in Ref. 



Ill ], and the dynamics of the early universe and the 



initial conditions for inflation in a model with radiation and a Chaplygin gas was considered 
in Ref. [3]. The main goal of the present work is to investigate the possible realization 
of a warm- Chaplygin inspired inflationary model, where the universe is filled with a self- 
interacting scalar field and a radiation field. We use astronomical data for constraining the 
parameters appearin g in this model. Specifically, the parameters are constrained from the 
WMAP observationsQ,Q. 

The outline of the paper is a follows. The next section presents a short review of the 
modified Friedmann equation and the warm- Chaplygin Inflationary phase. Section UTTl deals 
with the scalar and tensor perturbations. In Section IIVI it is presented a chaotic potential 
in the high dissipation approximation. Here, we give explicit expressions for scalar power 
spectrum and tensor-scalar ratio for our models. Finally, sectJV] summarizes our findings. 
We chose units so that c — h — 1. 



II. THE MODIFIED FRIEDMANN EQUATION AND THE WARM-CHAPLYGIN 
INFLATIONARY PHASE. 



We start by writing down the modified Friedmann equation, by using the FRW metric. 
In particular, we assume that the gravitational dynamics give rise to a modified Friedmann 
equation of the form 



H 2 = K [p ch + p 7 ] 



K 



1+/3 



+ P 7 



(3) 



where k = 87rG/3 = 8n/3mp (here m p represents the Planck mass), p^ = <fi 2 /2 + V(4>), 
V{4>) — V is the scalar potential and p 7 represents the radiation energy density. The 
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modification is realized from an extrapolation of Eg. ([2]), so that 



Pch 



[A + p£ +fr >] 



(1+/3)1 1+/3 



A + 



1+/9 



(4) 



where p m corresponds to the matter energy density [lOj. The generalized Chaplygin gas 
model may be viewed as a modification of gravity, as described in Ref. 3], and for chaotic 



inflation, in Ref. lOj. Different modifications of gravity have been proposed in the last few 



years, and there has been a lot of interest in the construction of early universe scenarios in 
higher-dimensional models motivated by string/M-theory 16|. It is well-known that these 
modifications can lead to important effects in the early universe. In the following we will 
take (3 = 1 for simplicity, which corresponds to the usual Chaplygin gas. 

The dynamics of the cosmological model in the warm- Chaplygin inflationary scenario is 
described by the equations 

4>+ 3H cP + V^-T 0, (5) 

and 

p 7 + AH Pl = T0 2 . (6) 

Here V is the dissipation coefficient and it is responsible of the decay of the scalar field into 
radiation during the inflationary era. V can be assumed to be a constant or a function of the 
scalar field 0, or the temperature T, or both [lj. Here, we will take T to be a function of 
only. In the near future we hope to study more realistic models in which V not only depends 
on but also on T, expression which could be derived from first principles via Quantum 



Field Theory [171, ll8(. On the other hand, T must satisfy r = f((f>) > by the Second Law 



of Thermodynamics. Dots mean derivatives with respect to time, and = dV{(f)/d(j). 

During the inflationary epoch the energy density associated to the scalar field is of the 
order of the potential, i.e. p$ ~ V, and dominates over the energy density associated to the 
radiation field, i.e. > p 1 . Assuming the set of slow-roll conditions, i.e. <fi 2 V(<f>), and 
<ti (3H + T)0 the Friedmann equation (j3J) reduces to 



and Eq. ([5]) becomes 



3#[l + r]0^ -V, 



where r is the rate defined as 



(7) 
(8) 



For the high (weak) dissipation regimen, we have r ^> 1 (r < 1). 

We also consider that during warm inflation the radiation production is quasi-stable, i.e. 
p 7 <C 4ifp 7 and p 7 <C Ftp 2 . From Eq.flB]) we obtained that the energy density of the radiation 
field becomes 

p 7 = ^ (io) 

which could be written as p 7 = aT^, where a is the Stefan-Boltzmann constant and T r is 
the temperature of the thermal bath. By using Eqs.(jE|), & an d fTTOT) we get 

r Vl 

rpi ,4> (~\~\\ 

Introducing the dimensionless slow-roll parameter e, we write 

" 1 <^ (12) 



if 2 6/c(l + r) (A + K 2 ) 3 / 2 
and the second slow-roll parameter 77 

# 1 



T/ , n 3 vv 2 

V.dxb + 



7 2(A + \/ 2 )_ 



(13) 



HH 3k (1 + r) vMTV 2 

where = d 2 V/d(j) 2 . 

Note that for r = (or T = 0), the parameters £ and rj given by Eqs.f ll2l) and ( jl~3l) 
reduced to typical expression corresponding to cool inflation where a the Chaplygin gas is 



considered 



ioj. 



It is possible to find a relation between the energy densities p 7 and p^ by using Eqs. ffTTj) 
and (fl2|) . so that 

r ( j4 + c«) ,,.,1 
^ = 2(l + r) " (14) 

Warm inflation takes place when the parameter e satisfying the inequality e < 1. This 
condition is analogue to the requirement that a > 0. The condition given above is rewritten 
in terms of the energy densities and p 7 , so that 

(A + p 2 ) 2(1 + r) 

> P 7 - I 15 ) 

P* r 

Also, inflation ends when the universe heats up at a time when e ~ 1, which implies 

^±f!L?(i±i) P7 , ( i6) 
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and the number of e-folds at the end of inflation is given by 

y/A + V 2 



N 



-3 K 



V, 



(l + r)d<f>'. 



(17) 



In the following, the subscripts * and / are used to denote to the epoch when the cos- 
mological scales exit the horizon and the end of inflation, respectively. 

III. PERTURBATIONS 



In this section we will study the scalar and tensor perturbations for our model. Note 
that in the case of scalar perturbations the scalar and the radiation fields are interacting. 
Therefore, isocurvature (or entropy) perturbations are generated besides of the adiabatic 
ones. This occurs because warm inflation can be considered as an inflationary model with 
two basics fields 3- m this context dissipative effects can produce a variety of spectral, 
ranging between red and blue j^, 19 L and thus producing the running blue to red spectral 
suggested by WMAP observations 14( . 

As argued in Refs. Q, Q , the density perturbation could be written as 8h = Scj). From 
Eqs.® and ©, this expression becomes 



4 



36 k" 
~25~ 



{A+V 2 ){l + r) 
V 2 



(18) 



In the case of high dissipation, the dissipation coefficient T is much greater that the rate 



expansion H , i.e. r = T/3H 3> 1 and following Taylor and Berera[22j, we can write 



2tt 2 



(19) 



where the wave-number k F is defined by kp = ^/TH/V = H y/3r > H, and corresponds to 
the freeze-out scale at which dissipation damps out to the thermally excited fluctuations. 
The freeze-out wave-number k F is defined at the point where the inequality V ss < FH, is 



satisfied 



22] 



From Eqs. (pj| and (USD ^ follows that 

18^ T r 



4 



25 vr 2 V\ 



[nrVA + V 2 ) 5/2 . 



(20) 



The scalar spectral index n s is given by n s 



din si 



dink ' 



r L , where the interval in wave number 



is related to the number of e-folds by the relation d\nk{4>) = —dN(<f)). From Eq. (j2"U|) . we 



get 

n s « 1 - [he- 27} - C] , 
where, the parameters e, and £, (for r ^> 1) are given by 
VVl 



1 



1 



and 



6k r (A + V 2 f/ 2 ' 
Vl 

C 



V 



3nry/A + V 2 

V s 



n 3 VV^ 



V 2(A + V 2 )_ 

v ,<t> r A 



(21) 



(22) 



(23) 



respectively. 

One of the interesting features of the observations from WMAP is that it hints at a 
significant running in the scalar spectral index dn s /dlnk = a s [jjj]. From Eq. (l2TT) we obtain 
that the running of the scalar spectral index becomes 



a. 



2e(A + V 2 
VVs 



[5e^-27^-Crf 



(24) 



In models with only scalar fluctuations the marginalized value for the derivative of the 

n 

spectral index is approximately —0.03 from WMAP five year data[l_4|. 

As it was mentioned in Ref. the generation of tensor perturbations during inflation 
would produce stimulated emission in the thermal background of gravitational wave. This 
process changes the power spectrum of the tensor modes by an extra temperature depen- 
dently factor given by coth(fc/2T). The corresponding spectrum becomes 

2 



A. 



16tt ( H 
2^ 



rat 



coth 



2T 



4k 



irmt 



(A+ V 2 ) 1/2 coth 



k 
2T 



where the spectral index n g , results to be given by n g 



In 



coth[fc/2r] 



(25) 

— 2e. Here, 



Og, 1WU!» .V, U , & XV^ ,Ug — ^ j, 

we have used that A 2 oc k ng coth[A;/2T] as described in Ref. 23] 

For r ^> 1 and from expressions ( 1201) and ( |25i) we may write the tensor-scalar ratio as 

V 2 



R{k) = 1 pjl 



T r r 5 /2 (A + V 2 )3/4 coth [ 2T 



(26) 



k=k t 



1 /2 

Here, 5h = 2 /5 and k* is referred to k = Ha, the value when the universe scale crosses 
the Hubble horizon during inflation. 

Combining WMAP observations Q with the Sloan Digital Sky Survey (SDSS) large 



scale structure surveys [24j, it is found an upper bound for R given by R{k* ~ 0.002 Mpc 



-v 



< 



0.28 (95%CL), where ~0.002 Mpc -1 corresponds to I = r k ~ 30, with the distance to 
the decoupling surface r = 14,400 Mpc. The SDSS measures galaxy distributions at red- 
shifts a ~ 0.1 and probes k in the range 0.016 h Mpc _1 < k <0.011 h Mpc -1 . The recent 
WMAP observations results give the values for the scalar curvature spectrum Pn{k*) = 
255 2 H {k if ) / A ~ 2.3 x 10~ 9 and the scalar-tensor ratio R(k*) ~ 0.055 < 0.2. We will make use 
of these values to set constrains on the parameters for our model. 



IV. CHAOTIC POTENTIAL IN THE HIGH DISSIPATION APPROACH 

Let us consider an inflaton scalar field <ft with a chaotic potential. We write for the chaotic 
potential as V = m 2 2 /2, where m is the mass of the scalar field. An estimation of this 



parameter is given for cool-Chaplygin inflation in Ref.[10j and for warm inflation in Ref. 22]. 
In the following, we develop models for constant and variable dissipation coefficient T, and 
we will restrict ourselves to the high dissipation regimen, i.e. r ^> 1. 



A. r = const. = To case. 

By using the chaotic potential, we find that from Eq.([H]) 

0=-^^ 0(t) = o e- m2t / r °, (27) 

i 

and during the inflationary scenario the scalar field decays due to dissipation into the radi- 
ation field. The Hubble parameter is given by 

/ 4 \ V4 

H(t) = ^[A + ^4e-^ 2 ^) . (28) 

Note that in the limit A — > the Hubble parameter coincide with Ref. 22|. The dissipation 
parameter r in this case is 

r / ™4 \ -i/4 

in 



where we observe that it increases when the cosmological time increases. 

By using Eq. (ll4p . we can relate the energy density of the radiation field to the energy 
density of the inflaton field to 

[Q fm 2 m p \ 

li \ 8IW pT^F 1 ' (29) 
8 



Note again that in the limit A — > 0, Eq. (l29|) coincides with that corresponding to the case 



where the Chaplygin gas is absent 



221 ]. i.e. p 7 oc 



1/2 



From Eq. fl20p . we obtain that the scalar power spectrum becomes 



9 T r 



An 2 m 2 V 



k 1 / 2 T (A + V 2 ) 1 ^ 



5/2 



(30) 



fc=fc» 



and from Eq. fT26l) the tensor-scalar ratio is given by 



R(k) 



rrir. 



3 V27T 



3 k 1 



1/2 \ 5/2 



m 2 V 



T T (A + V 2 yi 



coth — — 
2T 



By using the WMAP observations where Pk(K) ~ 2.3 x 10 9 , 
choosing the parameter T ~ T r , we obtained from Eqs. fl30|) and fl3T|) that 

5/4 



4 x 10" 5 T r 



^r 2 



coth(fc,/2T r ) 



(31) 
0.055, and 

(32) 



and 



A 



10 



-21 



coth 2 (fc,/2T r 



l/ 2 



(33) 



From Eqs. fl32|) and fl33|) . and since A > 0, we obtained a lower limit for m given by 

1/2 



m 2 > 10 6 T r 



p5 
1 



a/ coth[A;*/2T r 



(34) 



and also we find an upper limit for the chaotic potential when the universe scale crosses the 
Hubble horizon during inflation, given by 

,4 



K < 10 



-11 



ml 



(35) 



coth[A;,/2T r ] ' 

Now we consider the special case in which we fixe T ~ T r ~ 0.24 x 10 16 GeV and T ~ 
0.5 x 10 13 GeV. In this special case we obtained that the lower limit for the square mass of 
the scalar field, is given by m 2 > 4 x 10~ 15 m 2 , and for the chaotic potential we obtained an 
upper limit given by V* < 10~ 14 rrip. 

In Fig.(pQ) we have plotted the parameter (m/m p ) 2 versus the dissipation coefficient 
T /m p . In doing this, we have used Eqs. fl2~T|) . fl32l and fl33l) . Here, we have taken the 
values T ~ T r ~ 0.24 x 10 16 GeV and K = 0.002 Mpc" 1 . 



9 




FIG. 1: In this figure we plot the parameter (m/m p ) 2 as a function of the dissipation coefficient 
T /m p . Here, we have taken the values T ~ T r ~ 0.24 x 10 16 GeV and k* = 0.002 Mpc -1 . 



B. r oc 6 n case. 



In this case we consider a power-law dissipation coefficient given by 



T(cf>) = a n f 



(36) 



where n is a positive integer number and a n is a positive constant. The case n — 1 was 



studiec. 
in Ref. 



25 1 



in Ref. Q , the case n = 2 was developed in Ref. [20j , and the general case for any 



n 



From Eq.(fS1) we find that 



<j>{t) 



nm 



a. 



t 



l/n 



(37) 



where <fi(t = 0) 
becomes 



r(t) 



a 



The dissipation parameter, r, as function of the cosmological time 

2 \4/n1-V4 



3«V2 



m 



(38) 



10 



and in the high dissipation approach it is necessary to satisfy that 



a, 



3«V2 



> 



A + 



4 



,2 \ 4 /n' 
t 



a, 



1 1/4 



Using Eq. (1141) . we can relate the energy density of the radiation field to the energy density 
of the inflaton field as 



l-n/2 



7T V8 2«/ 2 a n y (A + p2)i/4- 
From Eq. fl20p . we obtain that the scalar power spectrum becomes 



(39) 



Pn(k) 



9VST r 



4^2 : m±± 



m 2 



2 n ' 2 H}/ 2 a n {A + V 2 f' A 



5/2 



(40) 



and from Eq. (l26l) the tensor-scalar ratio becomes given by 



R(k) 



m„ 



3V8tt V a 



3 K 



1/2 \ 5/2 



m 2 (2 1/) 



Tr (A + V2)i/8 coth ^2T 



(41) 



By using the WMAP observations where Pn(K) ~ 2.3 x 10 9 , R(K) = 0.055, and 
choosing the parameter T ~ T r , we obtained from Eqs. (l40p and (14 ip that 



3.7 x 10- 4 



m 
and 

A « 5.9 x 10~ n 
where a and h are given by 

2^ T r a 5 n /2 



avn 5 /*v* 4 

8.1 x 10~ 10 



2 

4+5n 



fi- 



fe 2 



- 1.7 x 10 10 K 2 



(42) 



(43) 



and 6 



m„3 3 / 2 k 5 / 4 , / jfe. 

coth 



2T r 



r- 2+5n 5/2 ^ 
VVT 2 4 an ±r 

respectively. 

Using that A > 0, then we find from Eqs. ( )42|) and ( 1431) a lower limit for the mass of the 
inflaton m, independent of the value V*, and it becomes 



m > 



,3-5n 



10 2 



2(4+5n) 



For instance, for n 

.511/9 



1 we obtained from Eq. 



(44) 



that m > 



10- 1 [T 2 < =l /^]i/9 [ mp coth(fc»/2T, 



r/J 1 / 3 ~ 7 x 10 7 m p . Here, we have used that 
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T r ~ 10~ 3 m p , K = 0.002MPC- 1 and a n=1 ~ 10~ 5 (see Ref.jgJ). Note that this lower limit for 
m decreases when the parameter a n= \ decreases. For n > 1 we noted that a similar situation 
occurs for the lower limit of m, i.e., this limit decreases when the parameter a n decreases. 
In FigfS] we have plotted the running spectral index a s versus the scalar spectrum index 
n s , for a dissipation coefficient T(0) = a n= \ <fi. In doing this, we have taken two different 
values for the parameter a n= \ and choosing the parameters T ~ T r ~ 0.24 x 10 16 GeV and 
fc* = 0.002Mpc~ 1 . Note that for a n= \ ~ 10 -5 — 10~ 7 and for a given n s ~ 0.95 the values 
of a s becomes in agrement with that registered by the WMAP observations. From Fig. [3] 
we note that the parameter A ~ 10~ 25 m p becomes smaller by one order of magnitude when 
it is compared with the case of Chaplygin cool- inflation [lOj, and smaller by two order of 
magnitude when it is compared with the case of tachyon-Chaplygin inflationary universe 



Another interesting situation corresponds to the case when V = a n=2 2 . Following a 



and A ~ 10 2& m^ in order to be in agrement with WMAP five year data[l_4j. Note that 
this value of A is similar to those found for the case n = 1. Here, again we have taken the 
parameters T ~ T r ~ 0.24 x 10 16 GeV, K = 0.002MPC" 1 and n s ~ 0.95. 

In FigJUwe show the plot of the adimensional parameter a n /m^~ n in terms of the param- 
eters n on the one hand, and m/m p in addition. In doing this, we have used Eqs. (121]) . fT24l) . 
( f42l) and ( f43l) . Here, we have used the WMAP five year data, and we have taken three differ- 
ent values for the parameters n and m/m p . Also, we have chosen T ~ T r ~ 0.24 x 10 16 GeV 
and = 0.002Mpc _1 , as above. 

V. CONCLUSIONS 

In this paper we have investigated the warm- Chaplygin inflationary scenario. In the slow- 
roll approximation we have found a general relationship between the radiation and scalar 
field energy densities. This has led us to a general criterium for warm-Chaplygin inflation 
to occur. 

Our specific models are described by a chaotic potential and we have consider different 
cases for the dissipation coefficient, T. In the first case, we took T = constant = T Q . Here, we 
have found solutions for the inflaton field and the Hubble parameter. The relation between 



model 11]. 



similar approach to those of the previous cases, we obtained that a n= 2 



10- 13 /m p -10- 5 



jm v 
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a 

s 




-0.03 -| 1 1 1 1 1 1 1 

0.90 0.92 0.94 0.96 0.98 




FIG. 2: Evolution of the running scalar spectral index a s versus the scalar spectrum index n s in 
the model r cx <f>, for two different values of the parameter a n= \. 

the radiation field and the inflaton field energy densities presents the same characteristic to 
that corresponding to the warm inflation case, except that it depends on the extra parameter 
A. For the case in which the dissipation coefficient V is taken to be a function of the scalar 
field, i.e. V oc n , it was possible to describe an appropriate warm inflationary universe 
model for n = 1 and n = 2. In these cases, we have obtained explicit expressions for the 
corresponding scalar spectrum and the runnin g of the scalar spectrum indices. 

Finally, by using the WMAP five year data|l4j and for the special case where T ~ T r ~ 
0.24 x 10 16 GeV, = 0.002Mpc _1 and n s ~ 0.95, we have found new constraints on the 
parameters a n , A and m. 
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A/m 

p 

6x10 26 H 1 1 1 1 1 1 1 r 




-| 1 1 1 1 1 1 1- 1 1 

0.90 0.92 0.94 0.96 0.98 1.00 




FIG. 3: The parameter A versus the scalar spectrum index n s in the model T oc 0, for two different 
values of the parameter a n= \. 
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